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In the present work we demonstrate that the spin columnar dimerized phase (conducting spin 
stripe liquid) is stable at doping x c i < x < x C 2- At x < x c i the system undergoes phase transition 
to the Neel state, and at x > x C 2 it is Normal Fermi liquid. At t/J = 3 the critical concentrations 
are x c \ ~ 0.09 and x C 2 ~ 0.36. To prove stability of the stripe phase and to calculate critical 
concentrations we employ two approaches. The first one is more technically involved, and it gives 
more accurate value of the critical concentration. The approach is based on the calculation of the 
magnon Green's function. Imaginary poles in the Green's function indicate transition to the Neel 
state. The second approach consists in direct comparison of ground state energies of the Neel state 
and the columnar dimerized state. Both approaches demonstrate stability of the spin stripe phase 
and give close values of the critical concentrations. 

I. INTRODUCTION 

It is widely believed that the 2D t — J model is relevant to the low energy physics of high-temperature super- 
conductors. This is why investigation of this model is of great interest both for theory and experiment. In spite 
of great efforts during more than a decade there is no full understanding of the phase diagram of the t — J model, 
however some facts are well establishecL, At zero doping the model is equivalent to the Heisenberg model on a square 
lattice which has long range Neel ordeiu. Doping by holes destroys the order. A simplified picture .of noninteracting 
holes leads to the Neel state instability with respect to spirals at arbitrary small but finite dopinga. However more 
sophisticated numerical calculations which take into account renormalization of the hale Green's function under the 
doping indicate that the Neel order is stable below some critical hole concentration x c ]P. In the_Neel phase (x < x c \), 
in all waves except s-wave, there is magnon mediated superconducting pairing between holes O. It is also clear that 
at very small hopping there is phase separation in the model because separation leads to reduction of the number of 
destroyed antiferromagnetic linkstj 

The purpose of the present work is to elucidate spin structure of the ground state at x > x c \ ■ The most important 
hint comes from experiment: indications of stripes in the high-T c materialsQ. Another important hint is a remarkable 
stability of the spin dimerized phase in the frustrated J\ — Ji model. The idf* of such state for this model was first 
formulated by Read and SachdevQ, and was then confirmed by further work&El The stability of such a configuration 
implies that the lattice symmetry is spontaneously broken and the ground state is four-fold degenerate. Such a route 
towards quantum disorder is known rigorously to take place in one dimension, where the Lieb-Schultz-Mattis (LSM) 
theorem guarantees that a gapped phase always breaks the translational symmetryE3. Some time ago Affleck suggested 
that the LSM theorem can be extended to higher dimensions, and the gapped states of quantum systems necessarily 
break the discrete symmetries of the latticeO. The example of the J\ — J2 model provides further support for this 
idea. 

There have been several attempts to consider the spin-dimerized phase in a doped Heisenberg antiferromagnet. For 
this purpose Affleck and Marstor£3 analyzed Hubbard-Heisenberg model in the weak-coupling regime, Grilli, Castellani 
and G. Kotlianl3 considered SU(N), N — > 00, t— J model, and very recently Vojta and SachdevEj considered Sp(2N), 
N — > 00, t — J model. These works indicated a stability of the spin-dimcrized phase in some region of parameters, 
providing a very important guiding line. However relevance of these results to "physical regime" of the t — J-model 
remained unclear. Stability of the spin-dimer order for the t — J model has been demonstrated in the paperEj. The 
only small parameter used in the analysis was hole concentration with respect to the half filling. In the present work 
we continue studies in the same direction applying various techniques. To be confident in the results we prove stability 
of the dimer phase by two independent methods: 1) Calculation of the magnon Green's function in the dimerized 
phase (Green Function Method), 2) Comparison of the ground state energies of the doped Neel state and the dimerized 
state (Direct Energy Method). The second approach is very simple physically and technically. The first approach is 
more technically involved, but it allows us to calculate the critical concentration more precisely. 

To incorporate some experimental data we consider t — t' — t" — J model defined by the following Hamiltonian 

(ij)a (ij'i)ff (i]2)<y (ij)<? 



1 



c itT is the creation operator of an electron with spin a (a =|, I) at site i of the two-dimensional square lattice. The 
(ij) represents nearest neighbor sites, (iji) - next nearest neighbor (diagonal), and (ij'2) represents next next nearest 
sites. The spin operator is Si = ^cl a a a pCif3 and the number density operator is rij = c lcr c iv The c\ a operators 
act in the Hilbert space with no double electron occupancy. Antiferromagnetic interactions > are arranged in 
a stripe pattern shown in Fig. 1: solid links correspond to Jy = J± = J(l + 5), and dashed links correspond to 
J ij =j = J(l-5). 
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FIG. 1. Stripe spin dimerization on square lattice. Solid links correspond to J± — J(l + 5), and dashed links correspond to 
j = J(l-5) 

For real cuprates the antiferromagnetic interaction is isotropic, 6 = 0. However for theoretical analysis of the 
magnon Green's function it is convenient tQ_consider nonzero 5 and later set 5 — ► 0. The antiferromagnetic exchange 
measured in two magnp» Raman scatteringc3 is J = 125meV. Calculation of the hopping matrix elements has been 
done by Andersen et aE3. They consider a two-plane situation, and the effective matrix elements are slightly different 
for symmetric and antisymmetric combinations of orbitals between planes. After averaging over these combinations 
we get: t = 386meV, t' = — 105meV, t" = 86meV. Below we set J = 1, in these units 

t = 3.1, £' = -0.8, t" = 0.7 (2) 

These values are confirmed by the analysisEl of photoemission (PES) data for charge transfer insulator S^CuG^C^. In 
further numerical estimates we will use the values (||) , but we will also consider "pure" t — J model which corresponds 
to t' = t" = 0. 

The rest of the paper is organized as follows. In Sec. II we describe magnon Green's function at zero doping, and 
remind ideas of the Brueckner technique used in the work. In Sec III the single-hole dispersion and wave function are 
considered. This section is important for both Green's Function Method and for Direct Energy Method. Sections IV 
and V contain the main results of the paper: In Sec. IV we demonstrate stabilization of the dimerization by doping 
using Green's Function Method, and in Sec. V we come to the same conclusion using Direct Energy Method. Sec. 
VI addresses the quantum phase transition from the dimerized liquid to the Normal Fermi liquid at high doping. In 
Sec VII we discuss shape of the Fermi surface and distribution of the photoemission intensity. Sec. VIII summarizes 
the work. Some technical details concerning so called "triple" diagrams are discussed in Appendix. 



II. ZERO DOPING 

At half filling ((n^) = 1) the Hamiltonian ([!]) is equivalent to a Heisenberg model which has already been studiedBEl: 
for 5 > S c ~ 0.303 the ground state is a quantum state with gapped spectrum, and for S < S c there is spontaneous 
Neel ordering with gapless spin waves. 

In order to study the stability of the dimer phase we first derive an effective Hamiltonian in terms of bosonic 
operators creating spin-wave triplets (magnons) t ia , a = x, y, z and fermionic operators creating holes aj., a =|, J, 
from the spin singlets shown in Fig. 1. This Hamiltonian consists of four parts: the spin-wave part H t , the hole part 
Hh, the spin- wave-hole interactioruiJt/i, and the hole- hole interaction H^h- Let us start from H t . Similar effective 
theories have been derived in Refs.t£l and we only present the result: 

H t =H 2 +H 3 + H 4 + Hu, (3) 

, H -i = R ( k l' k 2) e a/37 t k 1 «*k 2 A37 + h - C - 

1+2=3 
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Hi = 2_j T(ki — k 3 )((5 Q( 5<5 / 3 7 — 5 a i35 1 s)t\ !Lia t^pt\ s _ 3l t\ ii s- 

1+2=3+4 

We also introduce an infinite repulsion on each site, in order to enforce the kinematic constraint t\ a t\p = 0. 

H V = UJ2 AoApUoU* , U -> oo (4) 

The following definitions are used in (||): 

A k = J ± + B k , (5) 
Bk = j(cosk y — 0.5cos£: x ), 
T(k) = j (0.25 cos k x + 0.5 cos k y ), 
R(p>q) = 0.25j(smq x ~ sinp x ). 

Throughout the paper we work in the Brillouin zone of the dimcrizcd lattice. 

At zero doping = 1) H t is-an-exact mapping of the original Hamiltonian (Q). To analyze this case it is enough 
to apply the Brueckner techniqucEIrtl The result for the normal spin-wave Green's function reads: 

n n ^ w + A k (-w) 

Gjv(k,cj = = = =t- (6) 

+ A^-uj)}^ - A k (u;)} + Bl 

where 

i k H = A k + sf(k) + 4 1 )(k, w ), (7) 

B k {u) = S k + Sf(k). 

Normal and anomalous T,f self-energies are caused by the quartic interaction H4 and the most important contri- 
bution Tigl comes from the Brueckner diagrams as described in0. Strictly speaking there is also some contribution 
to the self-energy caused by the "triple" interaction H3. However this contribution is very small (see, e.g. RefB) and 
therefore we neglect it. 

Expansion of the self-energy in powers of uj near uj — gives quasiparticle residue and spin- wave spectrum 



(8) 



^ = Z k ^[A k (0)Y -B k . 



Expressions for effective Bogoliubov parameters u k and v k are given in0. The spin- wave gap A = co ko , ko = (0, 7r), 
obtained as a result of a selfconsistent solution of Dyson's equations is plotted in Fig. 2 (line at x = 0). The critical 
value of the explicit dimerization (point where the-gap vanishes) S c = 0.298 is in agreement with results of series 
expansionsB and quantum Monte Carlo simulationscJ. The validity of the Brueckner approximation is justified by the 
smallness of the gas parameter n t = ^2 a (t\ a tia)- At the critical point n t — 0.13. 
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FIG. 2. The magnon "gap" (centre of gravity of spectral function) as a function of explicit dimerization S for tj J = 3 and 
different values of hole concentration x, x = 1 — (m). 

III. SINGLE HOLE DISPERSION 

Consider now doping by holes. On the single dimer |s) the hole can exist in symmetric (bonding) and antisymmetric 
(antibonding) states. Corresponding fermionic operators b\ a , and a\ a u =],[ creating hole from the singlet \s)i arc 
defined as: 

6t| s ) = -L( c t )CT+c t CT) |o), (9) 



alls) = 



^=(4, CT -cUlo), 



where 1 and 2 numerate the dimer sites. 

We will see later that doping suppresses spin quantum fluctuations of the dimerized state so that n t = 
at x > x c does not exceed 0.07. This is why we can neglect these fluctuations and consider pure dimerized spin liquid. 
For comparison we can say that even for spin ladder influence of the quantum fluctuations on the hole dispersion is 
not strongc3cJ in spite of the fact that in this case n t ~ 0.3. 

In leading approximation the wave function of a hole with given quasimomentum is of the form 

=E e4krn& iTis>> (10) 



where 15*) is the spin dimerized state, index n numerates the dimers and N2 — N/2 is number of sites in the dimerized 
lattice. We remind that throughout the paper we work in the Brillouin zone of the dimerized lattice. Sometimes the 
results are transfered to the usual lattice, but then it is specially pointed. The dispersion corresponding to the state 
( |To| ) can be easily calculated considering all possible hoppings between the dimers. The result is 

ei(k) = (l|i?|l> =-t+(t + t') cosfcy + (^ + t")cosk x + t' cosk x cosk y + t" cos2k y . (11) 

There is also an additional t-independent constant in the dispersion 

e = 1.75J. (12) 
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The constant arises because a hole destroys one spin dimer (this gives 0.75 J) and four links 1/Amrij (see Hamiltonian 
(|l|)). In the present section we ignore en because it gives just a shift of the entire dispersion. However in Section V, 
calculating total energ y o f the system we will restore en. 

The wave function ( |l0| ) as well as the dispersion (|ll| ) are renormalized due to virtual admixture of antibonding 
states and magnon excitations. Now we calculate this admixture to show that it is small. In this calculation we follow 



the approach developed earlier for the doped spin laddered. The Hamiltonian ([j]) admixes following states to the wave 
function (llT 
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The state |2) is similar to |1) with replacement of bonding orbital to the antibonding one. The states |3) — 118) describe 
excited triplet (magnon) on the dimer n and a hole on the one of the closest dimers, for example n + x denotes the 
dimer on the right, n + x + y denotes the up-right dimer, etc. The brackets iya') denote that spins of the magnon 
and the hole are combined to the total spin 1/2 and z-projection 1/2: 1 1/2, 1/2). Calculation of matrix elements of 
the Hamiltonian is straightforward. The diagonal matrix elements are 
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The nonzero nondiagonal matrix elements are 
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Diagonalization of the Hamiltonian matrix can be performed numerically. This gives the quasiparticle dispersion ek 
and the quasiparticle residue which by definition is equal to the weight of the state |1) (see eq. ( |!o| ) in the exact 
wave function^.!. The minimum energy for different sets of parameters is shown in the fifth column of the Table I. 
We remind that we use units J = 1 . In the last column we show position of the minimum. For comparison we also 
display the minimum value of ei (see eq. ([il])), which always is at po = (ir, tt). 
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Let us denote the hole concentration by x = n/N , where n is number of holes, and N is number of sites. Hence on-site 
electron occupation number is (m) = 1 — x. Concentration of holes in terms of the dimerized lattice is two times 
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larger n/(0.5N) — 2x. At this stage we neglect interaction of holes between themselves, hence we consider them as 
an ideal Fermi gas, and the Fermi surface can be easily found from the condition 



21 W^' 



(16) 



where integration is performed inside the Brillouin zone of the dimerized lattice. The Fermi surface at x = 0.1 and 
hopping parameters given in (H) is shown in Fig. 3 by solid line. 
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FIG. 3. Fermi surface of the dimerized spin liquid at doping x — 0.1. Dot-dashed line corresponds to t/J = 3, t' ,t" = 0, solid 
and dashed lines correspond to the set of parameters from M). Solid line gives Fermi surface found with "exact" dispersion as 
it is described in Section III. Dashed line corresponds to the bare dispersion We put k x /2-K along the horizontal axis and 

ky/n along the vertical axis, so the picture corresponds to a quadrant of the Brillouin zone of the original lattice. 

We stress that in this figure we put k x /2Tr along the horizontal axis and ky/ir along the vertical axis. It means that 
the picture corresponds to a quadrant of the Brillouin zone of the original lattice. The Fermi surface corresponding to 
the bare dispersion (|ll|) is shown by the dashed line. The two curves are very close and this proves that the dispersion 
renormalization is small. The wave function renormalization is also not large and the quasiparticle residue is close to 
unity: it is = 0.83 at the bottom of the band, and at the Fermi surface Z^ h > = 0.80. So admixture of the states 
( |l3| ) to the bare state |1) is relatively small. In this admixture the states |3), |4), |7), and |8) clearly dominate. In a 
reasonable approximation the wave function can be written as 

^ « 0.9|1) - Q.22(|3) + |4) - |7> + |8». (17) 

Admixture of other components is even smaller. For comparison in Fig. 3 we show also by dot-dashed line the Fermi 
surface for "pure" t — J model (t/J = 3, t' = t" = 0, doping is the same, x = 0.1). We see that the additional 
hoppings influence substantially shape of the Fermi surface. 



IV. THE SPIN- WAVE-HOLE INTERACTION. STABILIZATION OF THE DIMER ORDER 

The magnon-holc interaction H t h can be easily calculated in the way similar to that for doped spin-ladder§lil. This 
interaction consists of two parts. The first one is interaction, pf a hole and a magnon positioned at different dimers. 
This is a relatively weak interaction which can be neglected£3'E2l. The second part, which gives the main effect, comes 
from the constraint that a hole and a magnon can not coexist at the same dimer: t[„ b]„ — t[„ a- = 0. To deal with 
this constraint we introduce, similarly to (0), an infinite repulsion 
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(18) 



The exact hole-magnon scattering amplitude caused by this interaction can be found via Bethe-Salpeter equation 
shown in Fig. 4a. 
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FIG. 4. (a) Bethe-Salpeter equation for hole-magnon scattering vertex V. Solid line corresponds to the hole and dashed line 
to the magnon. (b) Brueckner contribution to the magnon self energy. 

This scattering amplitude is similar to that for magnon-magnon scatteringS. Solution of the Bethe-Salpeter 
equation gives 



r(E,k) 



E 



-E-w q -ei(k 



(19) 



where E and k is total energy and total momentum of the incoming particles. It has been demonstrated in the 
previous section that the renormalized hole dispersion is close to the bare one t\ (k) . This is why we use in eq. ( fisj ) 
the bare dispersion (|TT|) . 

We remind that concentration of holes in terms of the dimerized lattice is 2x <C 1, and this is the gas parameter of 
the magnon-hole Brueckner approximation. In the previous section we have shown that the holes are concentrated in 
the pocket in the vicinity of po = (tt, tt). Therefore the magnon normal self-energy described by the diagram Fig. 4b 
is 



E^(k,c;) = 2xr[w + ci(p ),k + p ] 



(20) 



It is instructive to consider first the limit which allows an analytical solution: j <C J±, \[2txx -C 1. Bare magnon 
dispersion in this case is Wk ~ J± +j (cos k y — 0.5 cos k x ) and hence the integrals in (|9|,g0|) can be calculated analytically 
with logarithmic accuracy. This gives 



2V2Trx(t + j) 
ln(12.5/n) + iTT9(5tu)' 



where 



Slu = 
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t + j 



(21) 



(22) 



/i = max(|&j|, V2irx), 
and 9(5lo) is a step function. The magnon Green's function is 

1 



G(k,w) 



o - uj k - S^(k,w) 

For illustration the spectral function ImG(uj) at k = k = (0, n), t/j = 3 and different x is plotted in Fig. 5. 



(23) 
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FIG. 5. Magnon spectral density at k = ko = (0, tt) in the limit J± j, t/j — 3, and different hole concentrations. 

There are several conclusions from formula ( |2l| ) and Fig. 5: 1) doping pushes the spin- wave spectrum up, 2) the 
effect is increasing with hopping t, 3) finite width appears, 4) there is only a logarithmic dependence on the infrared 
cutoff. Let us stress the importance of the point (4). It means that the effect is practically independent of the long- 
range dynamics. Moreover, near the critical point (A = 0) the situation is even better: the spin-wave spectrum is 
linear and even the logarithmic divergence disappears. Thus in the 2D case there is separation of scales which justifies 
Brueckner approximation. If we tried to apply the described approach to the ID case (say a doped spin ladder) we 
would get into trouble: power infrared divergence appears in Brueckner diagram and hence there is no justification 
for gas approximation. Let us also comment on the point (3) (width). There is also a "triple" contribution to the 
magnon self-energy shown in Fig. 6a, b. 




a b c 

FIG. 6. (a) Normal "triple" contribution to the magnon self energy, (b) Anomalous "triple" contribution to the magnon 
self energy, (c) hole-magnon "triple" vertex. 



This is a long-range contribution, and it can be shown that it does not influence position of the critical point 
(A = 0) in linear in x approximation. This fact will be proved in Appendix and this is why here we neglect the 
"triple" diagrams. However note that these diagrams influence the width of the magnon spectral function. 

In the general case there are two contributions to the Brueckner self-energy: which is due to the magnon- 

magnon constraint, and Tr B ' r which is due to the magnon- hole constraint. To find the spia wave spectrum one has to 
solve selfconsistently Dyson's equation for Green's function (^), as it is described in RefEH. Results for the spin- wave 
"gap" A as a function of explicit dimerization S for different hole concentrations x and t/J = 3 are plotted in Fig. 
2. These curves are practically independent of the longer range hoppings t' and t" . Strictly speaking at x ^ the A 
is not a gap because of the large decay width. What we plot is the position of the center of gravity of the magnon 
spectral function. However at A — ► the width vanishes, and therefore the critical regime is uniquely defined. 

It is clear from Fig. 2 that at t/J = 3 and x > x c \ « 0.090 the "gap" remains finite even at S — 0. This is regime 
of spontaneous dimerization. Critical concentrations for other values of t/J are presented in Table II. 
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1 
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Xcl 


0.132 


0.106 


0.090 







Thus the doping stabilizes the dimerized phase. The larger the hopping 4, the stronger the effect of stabilization. 
(The same follows from eq. (f2l|).) This statement is true-only if t/J < 10. At t/J ~ 10 there is a crossover to 
quasiparticles with higher spin (hole-magnon bound states)c3 which indicates transition to the Nagaoka regime. The 
small parameter of the Brueckner approximation is concentration of holes in the dimerized lattice: 2x. Therefore a± 
t/J = 3 one should expect ~ 20% accuracy in calculation of x c \. Note that the value of x c i is close to that found inEl 
from the Neel state. 

Another important gas parameter is density of spin fluctuations nt- It also proved to be small: At the critical point, 
S — 0, x = x c = 0.09, the density is n t « 0.07. 



V. DIRECT COMPARISON OF GROUND STATE ENERGIES OF THE NEEL STATE AND THE SPIN 

DIMERIZED STATE 

In the previous section we have demonstrated quantum phase transition from the Neel state to the spin dimerized 
state at some critical hole concentration x c \. We calculated the magnon Green's function in the dimerized phase, 
and the transition point was identified as the point where the magnon gap vanishes. This is a rigorous approach 
which allows to determine x c i with relatively high precision (~ 20%), however it is rather involved technically. An 
alternative method is direct comparison of the ground state energies of these two states. This is a very simple method 
which does not require introduction of the explicit dimerization. So throughout this section 6 = 0. 

The energy p«r site for the undoped Neel state is —1.17 = —0.67 — 0.5, where the first term is Heisenberg energy 
(see e. g. Ref£j), and the second contribution comes from the — l/4n»nj term in the Hamiltonian (|]). We remind 
that in our units J=l. In this section we consider only "pure" 4 — J model,pi. c. 4' = t" = 0. Energy of a single 
hole injected into the Neel background is — 3.17£ 4- 2.834 27 + 1, see e. g. Ref.c3, where the last term comes from the 
—l/AniTij term in the Hamiltonian (H). Therefore in linear in x approximation energy of the doped Neel state is 

E Nee i/N = -1.17 + (-3.17* + 2.834 27 + l)x. (24) 

Energy of the undoped columnar dimerized state without account of quantum fluctuations is (—0.375 — 0.5) N, where 
—0.375 is Heisenberg energy and —0.5 comes from the — l/4riinj term in the Hamiltonian (|TJ). Quantum fluctuations 
push this energy down. Using perturbation theory one can find that in linear in n t (triplet density) approximation 
this shift is ~ — 0.5n t N. According to previous section at the critical point n t w 0.07 and therefore the shift is tiny. 
Energy of a single hole injected into the spin dimerized background has been found in section III, e = e + e m i,i, where 
eo is given by eq. (|l2|) and the values of e m in are presented in Table I. For 1 < 4 < 4 one can fit e m in as £ m in ~ —3.04. 
Altogether this gives following energy of the doped dimerized state (linear in x approximation) 

E m mer/N = -0.91 + (-34 + 1.75)x. (25) 

At x — the Neel state energy ( p4| ) is lower than that of the dimerized state (p5|). The critical concentration (the 
transition point to the dimerized state) is defined by the condition En eel — E^imer- This gives following values of 
the critical concentration: 



Table III. 



t/J 


1 


2 


3 


X C 1 


0.14 


0.11 


0.10 



The values of x& in Table III are somewhat overestimated. The matter is that the single hole energy for the Neel 
state is known with high precision, at the same time similar energy for the dimerized state has been found in Section 
III with trial wave function which contains only 18 components. The true energy is lower than-thc variational one. To 
estimate this uncertainty we refer to the doped spin ladder. There is a variational calculationES for this system which 
is similar to to the calculation in Sec. Ill, and there are also exact numerical simulationsEU. Comparison shows that 
the variational method underestimate t m in (see Table I) by 10-15%. Taking this as an estimate we should replace the 
term —34 in eq. ( p5| ) by —3.34. Then we come to following values of the critical concentration 



Table IV. 



t/J 


1 


2 


3 


X C 1 


0.12 


0.09 


0.08 
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Altogether results for x& presented in Tables II, and III, IV and derived by absolutely different methods are in 
remarkable agreement. This gives a very strong confirmation of the phase transition to the columnar spin dimerized 
state. Physical reason for stability of the spin dimerized state is especially evident after the energy considerations: 
this is the gain in the hole kinetic energy, it is "easier" to propagate in the dimerized background. 



VI. SPIN-DIMER ORDER PARAMETER AND TRANSITION TO THE NORMAL FERMI LIQUID AT 

HIGH DOPING. 



We have discussed the transition to the Neel state at hole concentration x < x c \. It is clear that at large x (x > x C 2) 
there is a 2nd order phase transition to the normal Fermi liquid. Let us define the spin-dimer order parameter as 

/0 =(S 2 S 3 )-(S 1 S 2 ), (26) 

where sites 1, 2, and 3 are shown in Fig. 1. For perfect dimerized state p = 3/4. There are two mechanisms 
for reduction of the order parameter. The first one is due to spin quantum fluctuations which approximately give 
p — > 3/4 — n t . The second mechanism is direct effect of doping. Naive estimate is p — > 3/4(1 — 2x), however the hole 
wave function (|l7j ) is slightly different from the bare one and because of this the coefficient in the naive formula is 
slightly renormalized. All together this gives 

P=- A [l-2x{l + &a 2 )]-n u (27) 

where a — 0.22 is the admixture coefficient in eq. (|l7j). This formula is derived in dilute gas approximation, i.e at 
2x,n t -C 1, however for an estimate we can extend it to large x. Setting n% ~ 0.05 we find from ( p7f ) that p vanishes 
at x C 2 ~ 0.36. We repeat that this is only an estimate because the approach assumes that 2x <C L 
The phase diagram of the t — J — S model at zero temperature is presented in Fig. 7 



1.0 
X 



• Normal Fermi liquid 



0.5 - 

x c2= 0.36 



Dimerized spin liquid 



5 C =0.298 




FIG. 7. Phase diagram of the t — J — S model (t/J = 3) in the plane doping (x) - explicit dimerization (5). 



Because of the mobile holes the dimerized spin liquid at x c \ < x < x C 2 is a conducting state. Stability of this state 
is a very robust effect because it is due to the high energy correlations (typical energy scale ~ 2t). There are also 
low energy effects with typical energy scale ~ 2tx which can lead to hole-hole pairing, charge stripes, etc. We do not 
consider these effects in the present work because they are secondary with respect to the main one: spin dimerization. 
However we would like to note that there is a simple mechanism for charge stripes induced by the spin dimers: Because 
of the anisotropic dispersion, see Fig. 3, the charge response is enhanced at some momentum p — {j> x , 0). The effect 
is very sensitive to additional hopping parameters t' and t" , they can further enhance or suppress the response. 
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VII. SHAPE OF THE FERMI SURFACE AND PES INTENSITY 



Shape of the Fermi surface for the dimerized state at the doping x — 0.1 and hopping matrix elements from (g) is 
shown in Fig. 3 by the solid line. In this figure we put k x /2ir along the horizontal axis and k y /ir along the vertical 
axis, where k is defined on the dimerized lattice. In terms of the Brillouin zone of the original lattice this is usual 
quadrant: < P x < n, < P y < tt. To distinguish between the dimerized lattice and the original one we denote by 
capital letters momenta corresponding to the original lattice: P x = k x /2, P y = k y . 

In a real sample there are domains with one stripe dimcrization and there are domains with stripes rotated by 90°. 
Therefore in an experiment one should see two superimposed Fermi surfaces. Corresponding picture for x = 0.15 is 
shown in Fig. 8. 

1.0 i ■ ■ ■ — i — i 1 ■ — i — . . . 1 



Py 

7t 




0.0 1 1 1 1 1 1 1 1 1 1 1 

0.0 0.5 D 1.0 

rx 

IT 

FIG. 8. Two superimposed Fermi surfaces of the dimerized spin liquid corresponding to different orientations of the stripes. 
The hole concentration is x = 0.15. The picture corresponds to a quadrant of the Brillouin zone of the original lattice < P x < n, 

< Py < 71" . 

The first impression is that it isA£ery much different from what is usually observed in angular resolved photoemission 
(PES) experiments, see e. g. Ref.E3 However let us calculate intensity of the photoemission. 
The photoeffect operator is of the form 



i = ^c u e iPr ", (28) 



where summation is performed over sites of the square lattice. Amplitude of the hole creation from the dimerized 
background (oriented as it is shown in Fig. 1) is equal to 

A = (s|6 T i|s) = cos(P x /2). (29) 

Here we have taken into account that according to considerations in Section III wave function of the hole in bonding 
state is practically unrenormalized. We stress once more that P is a quasimomentum corresponding to the original 
lattice. According to ( p9| ) intensity of PES spectra (/ oc A 2 ) is dropping quickly as P x is increasing. If one assumes 
that the t — J model originates from the single band Hubbard model, then the corrections of the order oit/U to the 
PES intensity can be calculated in a way suggested in RefJia. This gives 

i P cx ^cos y + ii C0S ^ + i cos y cos Py ) ■ ( ' 30 ' ) 

The Hubbard repulsion U is excluded from this formula using relation J — At 2 /U . According to (]3(]) the PES intensity 
is strongly asymmetric at the Fermi surface. For example at x = 0.15 the intensity at the right top corner of the Fermi 
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surface, P x = 0.667T, P y = tt (see Fig. 8) is 3.5 times smaller than that at the left top corner of the Fermi surface, 
P x = 0.347T, P y — 7T. In real cuprate the asymmetry must be even stronger. The reason for further enhancement of 
the asymmetry is structure factor of the Zhang- Rice, singlet. This is quite similar to the well understood situation in 
the charge transfer insulator Sr 2 CuC>2Cl2, see Ref.Ea. 

Thus the angle resolved photoemission measurements are sensitive only to the "inner" parts (the parts closest to 
the r point: P = (0,0)) of the Fermi surfaces shown in Fig. 8. Shape of this effective "Fermi surface" is very close 
to what is observed in numerous PES experiments. Another feature which agrees with experiment is width of the 
"quasiparticle" peak along (1,1) direction: it is always rather big because the peak arises as a superposition of two 
different peaks corresponding to two different Fermi surfaces. 



VIII. CONCLUSIONS 



In conclusion, using the dilute gas approximation we have analyzed the phase diagram of the t — J model and 
the stability of spin-dimerized phase. The main result of the work is phase diagram shown in Fig. 7. Without any 
explicit dimerization (5 — 0) the spin dimerizcd phase is stable at hole concentration x c \ < x < x C 2- At t/ J — 3 the 
critical concentrations are x c \ ~ 0.09, x C 2 ~ 0.36. At x < x c \ the system undergoes transition to the Neel state, and 
at x > x C 2 to the Normal Fermi liquid. 

To prove stability of the spin dimerized phase and to calculate critical concentrations we have used two independent 
approaches. The first one is based on the calculation of the magnon Green's function. The second approach consists 
in direct comparison of ground state energies of the Neel state and the dimerized state. Both approaches demonstrate 
stability of the spin dimerized phase and give very close values of the critical concentrations. 
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IX. APPENDIX: "TRIPLE" DIAGRAMS 



Purpose of the present section is to demonstrate that the "triple" diagrams shown in Fig. 6a, b do not influence 
position of the critical point x c \ found in Section IV. The "triple" vertex is shown in Fig. 6c. In the present section 
we completely neglect small renormalization of the hole wave function considered in Section III. Therefore, the initial 
state in Fig. 6c is given by ( ^0| ) and the final state is 

\f) = bit\ a \s) = ^xy pr ° 6 -£ eiqrm *-i s >- (31) 

n m 

Kinematic structure of the vertex is obvious 

r p , q = -iAil a [bla a b p+(1 ], (32) 

where a a , a = 1,2,3 is vector of Pauli matrices, and bp is the hole wave function in spinor representation. Direct 
calculation of the matrix element (f\H\l) and comparison with (p3) gives following value of the coupling constant 



A = l 
2 



(t + 2t cos p y ) sin p x + — sin q a 



(33) 



The normal "triple" magnon self-energy E3„(q, ui) is shown in Fig. 6a. Since we are interested in the critical point x c i 
it is enough to calculate the self energy at zero frequency and at the momentum where the spin-wave gap vanishes: 
q = qo = (0, 7r). Straightforward calculation of the loop at t' = t" = gives 

£ 3 n(qo,0) = -tx. (34) 

Unfortunately analytical calculation at the values of t' and t" from (0) is impossible because at these values the 
curvature of the hole dispersion along y-direction vanishes. However numerical calculation shows that at x = 0.1 — Q. 15 
and t', t" from (0) the "triple" self-energy is by a factor 1.5 smaller than one given by (pi[). Comparing the "triple" 
self-energy with the Brueckner one ( pfj| ) we find that the "triple" self-energy is by a factor 7 smaller. This is already 
enough to neglect the "triple" contribution. However we would like to demonstrate that suppression of the "triple" 
contribution is even stronger: its influence on the point of phase transition is exact(!) zero in linear in x approximation. 
This interesting fact is related to anomalous "triple" self energy shown in Fig. 6b. Simple consideration based on the 
structure of the vertex (^) shows that following exact relation takes place 

£ 3 a(qo,0) =-E 3n (q ,0). (35) 
According to (^) the magnon spectrum found without "triple" diagram is of the form 

w q cx Vi 2 - B 2 , (36) 

where A arises from the normal terms in the effective Hamiltonian and B arises from the anomalous terms. At the 
critical point A = — B and hence the excitation energy vanishes. With account of "triple" diagrams the relation ( |3^ ) 
should be rewritten as 



w q « - (/.?-• E Ja ) 2 . C-571 

It is clear that because of ( |35| ) the dispersion ( |37| ) vanishes exactly at the same point as (|36|). This proves our 
statement that the "triple" self energy does not influence position of the transition point. However away from the 
transition point when the spin- wave gap is increasing the "triple" self. energy is getting more important. Note that 
similar situation takes place with "triple" diagrams considered in Rcf.0 for J\ — J2 model. 
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